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Abstract. An approach to calculate radiative corrections to unpolarized cross section of semi-inclusive 
electroproduction is developed. An explicit formulae for the lowest order QED radiative correction are 
presented. Detailed numerical analysis is performed for the kinematics of experiments at the fixed targets. 
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1 Introduction 



Semi-inclusive processes of hadron electroproduction have 
been recognized long ago Q as an important tool for test- 
ing QCD predictions of nucleon structure because them 
admit to get information about quark distributions in 
the nucleon for each flavour separately. Precise analysis 
of the hadron structure functions extracted from the ex- 
perimental data requires, however, an iterative procedure 
involving radiative correction (RC) of these data. RC to 
cross section of semi-inclusive processes can be calculated 
on basis of Bardin and Shumeiko covariant approach of- 
fered in ref . [|| originally for elastic scattering. In ref. (3) the 
method was developed on semi-inclusive processes, where 
analytical formulae for the lowest order RC for coincident 
processes in electroproduction were found. In this paper 
we present analogous formulae, but in contrary to results 
of the ref.|| we do not assume integration over hadronic 
kinematical variables p\ and 4> h . It allows to calculate the 
model-independent RC relying only on the common rep- 
resentation for the hadronic tensor. 

Last years the cross sections of the hadron electropro- 
duction on fixed targets were measured as functions of 
azimuthal angles and transversal momentum of registered 
particles (see j| and references therein). However this in- 
formation is not sufficient to extract all structure func- 
tions involved in the hadronic tensor in wide kinematical 
region required for RC calculation. Therefore we have to 
use some model for the semi-inclusive structure functions. 
Such a model should give a good initial approximation for 
the iterative procedure, and we hope that the model for 
structure functions, which can be constructed on the basis 
of results of Mulders and Tangerman || , is an appropriate 
one. It should be noted that unpolarized structure func- 
tions were considered in [^) along with the spin dependent 
distributions. However, in this paper we restrict ourselves 
only by calculation of RC to unpolarized cross section, and 
the consideration of RC to observable quantities in polar- 
ization experiments on hadron electroproduction will be 
a subject of a separate publication. Notice that for az- 



imuthal effects we use additionally a model given in ref. q{ 
(see also 0). 

In the Section || we shortly describe the kinematics of 
the hadron electroproduction process with and without ra- 
diation of additional photon. The hadronic tensor and the 
model for structure functions are discussed in Section |^. 
The analytical formulae for Born cross section and for RC 
of the lowest order are given in Sections || and |]. In this 
paper the only ultrarelativistic approximation is made: 
electron mass is considered to be small. We note that fi- 
nal analytical formulae are written in the form similar to 
used in FORTRAN code POLRAD 2.0 §. The Section | 
is devoted to numerical analysis, performed on the basis 
of new code HAPRAD specially developed by us for this 
purpose. Most cumbersome formulae are gathered in the 
Appendix. 



2 The kinematics 

The cross section of hadron h electroproduction 

e(fci) + N(p) — » e'(fc 2 ) + h{ Ph ) + X{p x ) (1) 

depends on five kinematical variables which can be chosen 
as 

x, y; z, t, (j) h , (2) 

where x and y are usual scaling variables, z and t are 
defined via hadron momentum 

t={q~Ph) 2 , z=p h p/pq, q=ki-k 2 , (3) 

<ph is an angle between planes (ki,k2) and (q, ph) in the 
rest frame (p = (M, 0)). Also the following invariants will 
be used 



S = 2kjp, X = 2k 2P = (1 - y)S, Q 2 



xyS, 



W 2 = S x -Q 2 + M 2 , S X = S-X, S P = S + X, 



Ml 



--p 2 x = (l- z)S x + M 2 + t, Vi, 2 = 2k li2Ph . (4) 
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When the radiative process 

e(fci) + N(p) — » e'(fc 2 ) + 7 (fc) + h{p h ) + X(p x ) (5) 

is considered, the three additional independent variables 
have to be introduced 



R = 2kp, t = qk/kp, 4> k , 



(6) 



4>k is the rest frame angle between planes (ki,k2) and 
(q, k). Also we introduce the quantity fj, = kph/kp and 
the following invariants 

Q 2 = (q - kf = Q 2 + Br, 
W 2 = {p + q-kf = W 2 - R(l + r), 
t=(q-k- p h ) 2 = t + R{r- fi), 

M 2 =pl = M 2 + R(l + T-v). (7) 

The phase space of the three final particles is param- 
eterized as 

d 3 k 2 d 3 p h d 3 k S x dzdtd<j) h RdRdrd4> k 

nb x dxdy — - — j= j= — . (8) 



k 2 Q Pho k 



2JX 



Q 



Instead of t— dependence we will also consider the cross 
sec tion as a function of transversal momentum p t defined 
in 



We are interesting in explicit dependence on angles 4>h 
and <pk- So it is useful to take the some scalar products 
with ph in the form 



Vl 2 = hnPh 



6cos( 



(9) 



/iR = kph — R{a k + b k (cos (ph cos <j>k + sin <p h sin </>&)). 



Also we will use a = a 2 ± a 1 . The expli cit expressions for 
coefficients are given in Appendix (A. 7). 



3 The hadronic tensor and the Born 
approximation 

The cross section of the electroproduction process can be 
obtained in terms of convolution of leptonic and hadronic 
tensors. There are two leptonic tensors: with and with- 
out additional radiated photon. The Born leptonic tensor 
(without a photon) is standard, and radiative one is cum- 
bersome. The explicit expressions for the leptonic tensors 
and formulae for RC in terms of them can be found in 

IS 

The hadronic tensor without the T- and P-odd terms 
can be presented in the form [[ll[|l2|| 



(10) 

where 



~nu a V , <r<f ~ M . q pq ^ M . t Phq 

V =9 +-qT' P=P + ~qT' Ph=Ph + —Q2—' 



(11) 



and all of the SF depend on four kinematical invariants 
(for example, Q 2 , W 2 , t, z) The model for structure func- 
tions can be constructed on basis of results of the pa- 
per Mulders and Tangerman . Keeping only the leading 
twist contribution we have for structure functions 



H 2 



Pi 



f q {x)D q g, 



M 2 E 2 



ME h 



Y,e 2 q f q (x)D q g, 



and 



(12) 



(13) 



Q = Gi = bcxp(-bp t ), 

where b = R 2 / z 2 is a slope parameter and R is a param- 
eter of the model. 

The Born cross section has the following dependence 
on 4> h : 



dan 



N 



Co 



dxdydzdpfdtfth Q 1 



(A + cos (f> h A c + cos2<f> h A 0C ), 

— (14) 
where N = a 2 yS x / \/Aq- The coefficients A do not depend 

on <f>h more and they have the form 



A 



2Q 2 Hi 



n t- {SX - M 2 Q 2 )H 2 
+ (4aV + 2b 2 - M%Q 2 )H 3 
+(2Xa 1 + 2Sa 2 - zS x Q 2 )H 4 
A c = 2b(2a + H 3 + 5 P H 4 ), 
A cc = 2b 2 H 3 . 



(15) 



When integrated over the kinematical variables 4>h an d 
Pt the cross section (|l4|) coincides with the well known 
formula for semi-inclusive cross section calculated within 
QPM 



da 



2-na i 



xyz dxdydz yQ' 2 

Unfortunately, the p\ distribution 
1 da 



{y 2 + 2~2y)Y^e 2 q f q {x)D q . (16) 



J xyz 



p\ 



g, 



(17) 



calculated with the exponential slope (13) does not fit 
experimental data with sufficient \ 2 ■ So the more compli- 
cated model |j| with power dependence on p 2 seems to be 
more adequate. Another possibility to come to an agree- 
ment with the data consists in replacement of the Gaus- 
sian factor ( |l3| ) by the fit of experimental p\ -distribution 
taken in the form B 



g = g-i = 



i 



bz +p 2 



c+dz 



(18) 
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X 


y 


Q 2 

GeV 2 


2 


SIRAD 


HAPRAD 




without cuts 


with cuts 


Qi 


Si 


Q2 


Qi 


Qi 


Q2 


0.038 


0.677 


1.33 


0.25 


1.029 


1.033 


1.024 


0.982 


1.041 


1.025 


0.985 


0.062 


0.567 


1.82 


0.35 


0.996 


0.989 


0.989 


0.947 


0.989 


0.980 


0.951 


0.092 


0.529 


2.52 


0.45 


0.970 


0.961 


0.961 


0.934 


0.961 


0.956 


0.936 


0.131 


0.499 


3.38 


0.55 


0.945 


0.936 


0.933 


0.912 


0.934 


0.931 


0.906 


0.198 


0.476 


4.88 


0.65 


0.918 


0.902 


0.902 


0.889 


0.897 


0.897 


0.881 



Table 1. The results for RC factors to three dimensional semi-inclusive cross section obtained using FORTRAN codes SIRAD 
and HAPRAD (see text for further explanations). Kinematical points are taken from the Table 1 of ref . |l7j] . 



4 The radiative correction of the lowest order 

The cross section that takes into account radiative effects 
can be written as 



o bs = a e 5 ' nI (\ + 6 VR + S vac ) + a F . 



(19) 



Here the corrections <5; n / and S vac come from radiation 
of soft photons |Q and effects of vacuum polarization^. 
The correction Svr is infrared free sum of factorized parts 
of real and virtual photon radiation. These quantities are 
given by the following expressions 



ovr 

Sinf 



ir\2 r 

-a - 

7T 

glept , g. 

w i>ar 1 w ' 



1 



X' 



-2- - In — 



S'X'-Q 2 p\ 
Li 2 - 



l)ln 

hadr 
vac i 



2 S' S'X 
S'X' 



TT 

~6 



(20) 



where S' = X + Q 2 -V 2 , X' = S-Q 2 -^i, l m = lnQ 2 /m 2 
and Li 2 is Spence function or dilogarithm. 

The contribution of radiative tail has the standard 
form |J{n| 



OF 



TV f f 

j d<p k / dr^^%(r,^ fc )x 



dRR j - 



i i=i 
Hi 



(Q 2 + Rt) 2 



' Q 4 



■(21) 



Here 2M 2 t, 



max . mm 



— S x ±y/\Q and R 



(M 2 -(M + 



m n ) 2 )/(l + T — [i), 5j = l for j~l and 5j=0 otherwise. The 
explicit formulae for functions 8(T,(pk) can be found in 
Appendix. The structure functions have to be calcu- 
lated for Born kinematics, but Jii is calculated in terms 
of tilde variables (]?]) . 



5 Numerical analysis 

In this section we give numerical results for RC to semi- 
inclusive unpolarized cross section. For all cases RC factor 

1 There are explicit formulae for leptonic contribution to vac- 
uum polarization effect (see jl5| for example) and parameteri- 
zation of hadronic one [[16[ . 



is defined as a ratio of observed to born cross sections. Also 
we will speak about relative RC (or simply RC), which 
is difference between observed and born cross sections or 
asymmetries divided on Born ones. 

For definiteness we choose the kinematics of experi- 
ment HERMES at DESY. First, HERMES is the modern 
current experiment with rich possibilities for studies in 
semi- inclusive physics (see JlTj). Second, HERMES is the 
experiment with electron beam, so the relatively large RC 
in respect to muon DIS experiment are anticipated. 

The dependence of RC on x, y and z is widely discussed 
in ref. || , where quark-parton model was assumed. Similar 
results can be obtained using the formulae of this paper 
after integration over p t and 9h- In this paper we concen- 
trate on comparison of the codes constructed on basis of 
the two sets of formulae and on studying of the effects 
beyond the quark-parton model: azimuthal asymmetries 
and dependencies on transversal momentum p t . 



5.1 FORTRAN codes POLRAD 2.0 and HAPRAD 

The special FORTRAN code HAPRAD was developed to 
calculate the RC to five-dimensional cross section 
d 5 a/dxdydzdp 2 ddh- From the other hand there is the code 
POLRAD 2.0 with special patch SIRAD, which calculated 
RC to semi-inclusive three dimensional cross section ob- 
tained in QPM. 

In this section we show that numerical results for RC to 
da / dxdydz reproduced by these two codes coincide with a 
good accuracy. It can be seen from the Table [j]. In the ta- 
ble we represent RC to the cross section as it follows from 
the runs of the codes POLRAD 2.0 and HAPRAD. Since 
HAPRAD allows one to take into account the kinematical 
cuts and to use different models for p\ -slope, three fits for 
p\ -distribution and the cases with and without experimen- 
tal cuts were considered. The first fit for p\ -slope is defined 
in ( |i"3| ) , while the second and third ones are our fits of ex- 
perimental data Q using exponential (Q[ — Q\ at b = a/ z) 
and power ( |l8| ) functional forms. As the kinematical cuts 
on 4>h and pt of the measured hadron we took HERMES 
geometrical ones Jl8| . We can conclude from this analy- 
sis that neither important differences between SIRAD and 
HAPRAD results, if exponential model for p\ -distribution 
is used, nor dependence on slope parameter model and ap- 
plying of geometrical cuts are found. However RC takes 
a negative shift in the case of the model based on power 
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2 z 



Fig. 1. Radiative correction to the semi-inclusive cross section 
for kinematics of HERMES; \f§=7.19 GeV. Symbols from top 
to bottom correspond to the x=0.05, 0.45 and 0.7. The results 
for x=0.15 are skipped, because they practically coincide with 
ones for x=0.05. 



functional form (18). As it is discussed below RC depends 
on steepness of pfdistribution. It is a reason why models 
like 8{j%) (QPM, POLRAD 2.0) and © give larger RC. 
Within practical RC procedure in concrete measurement 
of do j dxdydz the model can be fixed only if the informa- 
tion about pj-distribution is considered additionally. 

Also two models for fragmentation function were con- 
sidered: simple parameterization of the pion data and 
modern model in the next-to-leading order QCD [ J21| . RC 
factor calculated using these models differs on several per 
cent. However this model dependence is less important be- 
cause it can be eliminated by applying an iteration proce- 
dure, where fit of extracted data is used for RC calculation 
in subsequent step. 



Fig. 2. Radiative correction to semi-inclusive cross section vs 
pt; VS=7.19 GeV, x=0.15,Cf=4 GeV 2 . Dashed and solid 
curves correspond to models ( [j^ j and |71| ) respectively. 



the experimental fit for Q is used (solid curves in the Fig- 
ure |2|) there are no such a rise of RC for high pi values. 

The crucial for QCD predictions |l| quantity <pl> is 
expressed in terms of a as 



<Pt>~- 



I dp 2 t P 2 t a 
jdp 2 t d 



(22) 



RC to this quantity can be expressed as 



S P t = 



I d Pt Pt Vobs I [dpjjjj_Obo 

J dp} a obs I 
J dp 2 t p\ d obs 



J dp} o born 



J dp} a bo 



I dpi Pt &born 



J dpi a obs 



Opt 

5 ' 



(23) 



where 6 is semi-inclusive RC factor discussed above (see 
Figure |l]). The correction S pt both for exponential and 



5.2 Cross section and <pl> 

Here we give numerical results for unpolarized cross sec- 
tion in kinematics of experiment HERMES . RC factor 
(5) to the semi-inclusive cross section integrated over 
and <fih as a function of x, y and z is presented in Figure 
Further we analyse the z and pt dependence of the cross 
section and azimuthal asymmetries. The dependencies of 
RC factor (<5 = cr obs /a born ) to the semi-inclusive cross sec- 
tion on hadronic variables z and pt are shown in Figure 
where sigma bar (tr) is meant the four-dimensional cross 
section do j dxdydzdpl . We note that the obtained large the experiments |19 
correction to the cross section vs pt is an analog of the 
similar results for vector meson electroproduction [[lOj . In 
our case the slope parameter depends on z (see Q13D), so 
we have important z-dependence of the effect. However, if 



power model (jl^jlj) is presented on Figure [_ 

Similar to the case of pi -distribution the radiative ef- 
fect is larger when exponential model ( |l3| ) is used. That is 
because of contribution of small pi when we integrate over 



the phase space of emitted photon (21). The steeper the 
slope of the p^-distribution the larger t 
RC. 



lis contribution to 



5.3 Azimuthal asymmetries 

The following azimuthal asymmetries are measurable in 



2n 



< COS (f) h > = a" 



t> h cos(p h a obs , 
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Fig. 3. Radiative correction to <Pt> defined in eq. (pffi) for 
HERMES kinematics, VS = 7.19 GeV, y=0.4. Curves from top 
to bottom corresponds to x=0.15, 0.05 and 0.45. 

< cos 2 <j) h > = a~ b \ J d<j) h cos 2<j) h a obs , 
o 

< sin <f> h > = a~ b \ J d<p h sin <j) h cr obs . (24) 
o 

In terms of these quantities the observed cross section ( [l9| ) 
can be written as 

Cobs = &obs (l+ < COS (j>h > COS <f>h + 

+ < cos 2(j> h > cos 2<t>h) + Uadd- (25) 

Here o~ a dd is originated from contribution of higher har- 
monics (sin^, sin 20ft, ...). There are no their contribu- 
tions at the Born level (see eq.|l4]), and o-^J n = 0. 

Azimuthal asymmetry < cos 4>h > is negative in con- 
sidered region. RC to the quantity can exceed 10% The 
result is done in FigJJ. 

Within the model ( |l2] ) <cos2^; l > is equal to zero at 
the Born level. So the asymmetry in this case is defined by 
RC only. Our estimation shows that this effect is of order 
1%. Relative RC to the asymmetry can be estimated using 
another model ||, where <cos 2</>/j>^ at the Born level. 
It is of order 10% in the region of applicability of the 
model. 

< sin (j)^ > is equal to zero at the Born level in any 
case. But there is nonzero contribution to it coming from 
RC. The numerical analysis shows that radiative correc- 
tions does not give visible contribution to it. For HERMES 
kinematical region values of < sin <fih > does not exceed 
0.01%. 

We should stress that our predictions for the values of 
the radiative effects display a strong model dependence. 



Fig. 4. Azimuthal asymmetry <cos 4>h> vs z for y =0.2 within 
HERMES kinematics; \f§=7.19 GeV. Dashed (solid) lines cor- 
respond to born(observed) asymmetries. Curves from top to bot- 
tom correspond to x=0.1, 0.45 and 0.05. 



Therefore any reliable method of radiative correction of 
experimental data has to be based on an iterative proce- 
dure, where all necessary fits for RC codes use the pro- 
cessing data as an input and are specified in every step of 
this procedure. This procedure can be readily developed 
on the basis of the code HAPRAD. 



6 Discussion and Conclusion 

In this paper the QED radiative correction to different 
observable quantities in the experiments on hadron elec- 
troproduction is analyzed. The explicit covariant formulae 
are given in Section and Appendix. 

New FORTRAN code HAPRAD is developed in order 
to p erform the numerical analysis. It was shown in Section 
|5.1| that the results for the RC to cross section integrated 
over p t and (f>h are in agreement with POLRAD 2.0 
Several models for structure functions and slope param- 
eter in respect to p% were applied for. It is found that 
the model based on power p^-slope model leads to smaller 
values for RC. 

Within the exponential model for Q (|l^) RC to <pl> 
can exceed 40%. However it essentially depends on the 
model of pi -distribution as well as on x and z. 

RC to azimuthal asymmetries is of order 10%. The 
asymmetry <sin</>/j> due to RC is found to be negligible 
but not equaling to zero exactly. 

FORTRAN code HAPRAD is available (aku@hep.by) 
for the calculation of RC to observable quantities in the 
experiments on the hadron electroproduction. 
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Appendix 

In this appendix we list the explicit form for functions Oij : 

e l . J {T,(j) k ) = e%+cos(t) h e c iJ 

where 

e° 12 = 4tF ir , 



99 



y 43 



2FzS x - F d a k S x r + F d a~r + F d r 2 zS x 
—Fi + a k S p — Fi + a + , 

-F d cos <j) k b k S x T - F 1+ cos <j) k b k S p - 2F 1+ b, 



6l 3 = -b k sin 4>k{F d S x T + F 1+ S p ) 



J 4?, 

Here 



(A.2) 



?-+sin<M?.+cos2cMi/ ! (A.l) 



F F 

Fi+ = 1 , F 2 + — F I — 2- 

z\ z 2 \ zi, 



+F 2 S p + 2F 2+ M 2 t 
-2F ir M 2 t + F ir S x , 

= 2{-2F d b 2 T-F d {a+f 

_ I. _ _ 



":V2 



(tec 



#33 



6>g 2 = 2(-2F d b 2 r - F d (a+) 2 T - F 1+a - a 
+2F 2 _ cos 4> k bb k + 2F 2 -a k a+ 
-F IR M 2 T-2F IR a k a-), 



= 4(~2F d ba + T - F 1+ ba + F 2 - cos (j> k 
+2F 2 -ba k - cos <f) k F IR b k a~), 



.,b k a^ 



9 A2 



= Ab(-F d br + F 2 - cos <f> k b k ), 

= 4b k sin (f> k {F 2 -a+ - F In ar), 

= -2F d a + S p T - Fi + a~S p + Fi + a + S x 
+2F 2 -a k S p + 2F 2 _a+ + 2F 2+ tzS x 
+F IR a k S x - 2F IR a- - 2F ir tzS x , 

9% 2 = 2(-2F d bS p T + F 1+ bS x + F 2 _ cos cj) k b k Sp 
+2F 2 _b + cos (j) k F IR b k S x ), 

9t 2 = 2b k sm0 k {F2-S P + F IR S X ), 
0° 13 = -2(F + F d r 2 ), 

6° 23 = 2FM 2 + F d M 2 r 2 - ^F d S x r - ^F 1+ S p , 

6»° 3 = 2FMl + F d M%r 2 + 2F d a k a- T 
-2F 1+ cos (j) k bb k - 2F 1+ a k a + , 



= 2(F d cos <j) k b k a t - F 1+ cos <p k b k a + - 2F 1+ ba k ), 



OH = -2F 1+ cos 4> k bb k , 

6 S 33 = 2b k s\ncp k {F d a-T - F 1+ a+), 



F 2 _ = 



m 2 m 2 
zk z? 



F d = 



F 

z\z 2 '' 



where F = 1/(27twXq) and 

Fi R = F 2+ — Q 2 F d . 



(A.3) 



(A.4) 



The quantities z\ >2 — kkx^jkp in terms of integration vari- 
ables can be expressed as 

zi = X Q 1 (Q 2 S P + t{SS x + 2M 2 Q 2 ) - 2Mcos0fcV^A), 
z 2 = X q 1 (Q 2 S p + t(XS x - 2M 2 Q 2 ) - 2Mcos(j) k ^/KX), 

(A.5) 

where 

A r = (T-T min )(T max -t), X = SXQ 2 -M 2 Q 4 -m 2 X Q . 

(A.6) 

The scalar products of ph (M are expressed via coeffi- 
cients a 1 , a 2 , b, a k and b k : 

2Ma l = SE h - (SS X + 2M 2 Q 2 ) Pl X Q 1/ \ 
2Ma 2 = XE h - {XS X - 2M 2 Q 2 ) Pi X Q 1/2 , 
b 



-Pt\Jx/X Q , 



2Ma k = E h - pi(S x - 2M 2 t)X c 1/2 



b k = —Mpf\J X T / Xq. 

The quantities Eh, Pi and Pj are invariants 

zS x 



(A.l) 



E h = zv 
Pi\fX~ 



2M' 

r2 i /->2 



M 



= t - Mi + Q 2 + 2vE h , 
p 2 = El-p\-Ml. 



(A. 



In the rest frame they take sense of energy, longitudinal 
and transversal momenta of final hadron. 
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